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Abstract. Results of the study of the symbiotic binary Z And during its recent active phase 2000-2010 when it experienced 
a series of six optical outbursts are presented. High-resolution spectra obtained during the first and fourth outburst, which was 
the strongest one, have been analyzed. These data are compared with results of theoretical computations. The comparison 
provides information about the behavior of the system during the entire active phase rather than during an individual outburst. 
In particular it was found fundamental difference between the first outburst which opened the active phase and the recurrent 
outbursts — namely, the presence of bipolar collimated optical outflow during some of the recurrent outbursts. A scenario that 
can explain all the spectroscopic phenomena observed during this active phase as well as previous active phases of Z And is 
proposed. The possibility to use this scenario for explanation of the line spectrum of other classical symbiotic stars during 
their active phases is motivated. 
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1. INTRODUCTION 

Symbiotic stars are interpreted as interacting binaries consisting of a cool visual primary and a hot compact secondary 
component accreting matter from the atmosphere of its companion. Their spectral variability is determined from the 
orbital motion and the outburst events of the hot component. These events are often accompanied by intensive loss 
of mass in the form of optically thick shells, stellar wind outflow and bipolar collimated jets. The stellar wind is 
possible to be bipolar and highly collimated too. The ejected material forms complicated flow structure in the space 
between the components, whose elements determine the observed properties of the star. The interacting binary Z And 
is considered as a prototype of the classical symbiotic stars. It consists of a normal cool giant of spectral type M4.5 [1], 
a hot compact component with a temperature of about 1.5 x 10 s K [2] and an extended circumbinary nebula formed 
by the winds of the stellar components and partly photoionized by the compact object. 

Z And has undergone several active phases consisting of repeated optical brightenings with amplitudes up to 2-3 
mag and characterized by intensive loss of mass [3, 4, 5, 6, 7, 2, 8, 9, 10]. The last active phase of Z And began at 
the end of August 2000 [11] and continued up to that time including six optical brightenings. The maxima of the light 
during the active phase were in December 2000 (V ~ 8™8), November 2002 (V ~ 9 m 8), September 2004 (V ~ 9!"1), 
July 2006 (V ~ 8™6), January 2008 (V ~ 9™5) and January 2010 (V ~ 8!M) [9] (we also used data of the American 
Association of Variable Star Observers, AAVSO, and via private communication from A. Skopal). High resolution 
optical data indicating highly different physical conditions (velocity, density, temperature) in the line emission regions 
were obtained during this active phase [7, 2, 12, 8, 9, 10]. 

The current study is based on observations acquired during the first and fourth outbursts. The main features of the 
flow structure in the system during the active phase indicated by our data can be summarized as follows: 

1 . We suppose that an accretion disk is present in the system. Only in this way we can explain the presence of the 
He II A 4686 line in the spectrum and its behavior during the first outburst [8]. 

2. Stellar wind was observed in the system during all the outbursts. Two wind components with different velocity 
regimes were observed simultaneously — wind propagating at a moderate velocity, 0-100 kms _1 [2, 7, 8], and 
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FIGURE 1. Left panel: Schematic model of the region around the hot component during the first outburst. Right panel: The same, 
but in the plane perpendicular to the orbital plane where the emission regions are shown. (From spectroscopy presented in Tomov 
etal. [8].) 



wind with a high velocity, in excess of 500 km s- 1 [7, 8]. 
3. Bipolar collimated outflow appeared during the fourth outburst. In the beginning of July 2006 the H a line had an 
absorption component shifted by — 1400kms _1 from its center. It went into emission and later, in July-December 
2006, the spectrum contained additional emission components on either side of the central peak corresponding 
to velocities of 1200-1500 km s [13, 12, 9]. The line had the same features [12, 9]. Both the emission and 
absorption high-velocity components are assumed to be formed in a bipolar collimated outflow [12]. 

The main aim of the current study is to suggest a model for the flow structure in the system to explain all of these 
features. We believe that it will be helpful in explaining not only the data acquired during this active phase, but also all 
spectroscopic data taken during earlier phases, those following 1939, 1960 and 1984 and possibly the activity of other 
classical symbiotic stars. 

2. THE MODEL OF THE FLOW STRUCTURE IN Z AND 

According to the theoretical models [14, 15, 6] an accretion disk is formed around the compact object in the Z And 
system when the wind of its giant has quiescent parameters. The radius of the disk is about 5QRq for a wind velocity 
of 20-25 km s -1 and thus the outer part of the disk is optically thin. The mass of the disk is estimated as product of 
one quarter of the mass-loss rate of the giant ^2 x 10~ 7 Mq yr _1 [16], and the typical time interval between the active 
phases — about 10 years. Based on the size and the mass of the disk of ^5 x 10~ 7 Mq we can assume that its inner 
region is optically thick in the quiescent state of the system. 

To change the system from quiescent to active state a sufficiently large increase of the accretion rate is needed. 
Accretion of a considerable fraction of the disk's mass is required in order for an outburst to develop even in the 
combined model where the increased nuclear burning rate is taken into account [6]. Maximal increase of the accretion 
rate is possible in the framework of the mechanism proposed by Bisikalo et al. [14] and Mitsumoto et al. [15]. 
According to this mechanism, even a small increase of the velocity of the wind of the donor is sufficient to change the 
accretion regime. During the transition from disk accretion to accretion from the flow, the disk is partially disrupted 
and the increased velocity of the wind causes falling of the material of the disk onto the accretor's surface. However, 
even in this case a considerable amount of mass (up to 50-80 per cent) stays in the disk. Then massive accretion disk 
exists in the system during the active phase too. 

During the outburst the high velocity wind of the hot compact component collides with the accretion disk, its 
velocity decreases near the orbital plane but does not change at higher stellar latitude. A consequence of this is that 
two different kinds of stellar wind are observed — P Cyg wind with a low velocity and an optically thin high velocity 
wind. As a result of the collision the density of the wind close to the orbital plane increases and the level of the observed 
photosphere locates at a larger distance from the star. Thus, an optically thick disk-like shell forms, which occults the 
compact object. Since its effective temperature is lower it is responsible for the continuum energy redistribution and 



the increase of the optical flux of the star. The interaction of the wind with the disk is equivalent to interaction of two 
winds and thus can form shocked region of collisional ionization with temperature of up to 10 6 K (Fig. 1) [17, 6]. 

During the active phase the wind from the hot component "strips" the accretion disk carrying some part of its 
material in the circumbinary envelope. After each outburst some part of this material does not leave the potential well 
and after the cessation of the wind begins to accrete again. Because of the initial amount of angular momentum a 
disk-like envelope extending to larger distance from the orbital plane than the accretion disk itself forms (Fig. 2). The 
existence of a centrifugal barrier creates hollow cones around the axis of rotation with small opening angle [18, 19]. 
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FIGURE 2. The same like in Fig. 1 but during recurrent strong outburst. 

The disk-like envelope does not exist during the first outburst of the active phase. During the recurrent outbursts 
the extended disk-like envelope can collimate the wind which escapes the compact object predominantly via the two 
hollow cones (Fig. 2). The bipolar outflow can be observed as high-velocity satellite components situated on either 
side of the mean peak of the line. Their presence in the spectrum depends on the density of the disk-like envelope (i.e. 
on the system's activity during its previous phases) and on the wind intensity. They can appear only if the density of 
the disk-like envelope is sufficient to provide collimation and the mass-loss rate of the compact component is high 
enough. According to our model these spectral features can be observed during outbursts accompanied by loss of mass 
at high rate and preceded by similar strong outbursts. 

According to the current theory the existence of bipolar collimated jets is supposed to be due to presence in the 
system of magnetized disk which transforms the potential energy of the accreting material into kinetic energy of 
the outflowing jets. This means that accretion must be realized to provide existence of the jets. In the spectrum of 
Z And during its 2006 outburst, however, different indications of stellar wind presented along with satellite emission 
components. If we believe that "traditional" jets were observed, it is necessary to explain accretion onto the star at the 
time when a strong mass outflow from its surface was observed. On the other hand, the proposed model of collimated 
stellar wind makes it possible to explain all features of the observed spectrum without any important contradictions. 

Another question related to the nature of the collimated outflow is that about the orbit inclination. The velocity 
of the radiatively accelerated wind cannot exceed 3000 kms -1 . The highest observed outflow velocity during the 



2006 outburst was 1500 kms . A preliminary analysis of our spectra obtained during late December 2009 showed 
velocities of 1700-1800 km s _1 . Taking into account these velocities we obtain an upper limit of the orbit inclination 
of about 55°. That is why the view about the nature of the outflow can shed some light on the question about the orbit 
inclination. 

We suppose that the model suggested by us can be used to explain the behavior of the line spectrum of other classical 
symbiotic stars during their active phases too since the classical symbiotic stars have the next general characteristics: 

1. The mass transfer in the majority of the classical symbiotic stars is realized by means of the stellar wind of the 
giant. The theoretical computations of Bisikalo et al. [14] and Mitsumoto et al. [15] show that for systems with 
parameters close to those of Z And an accretion disk from wind accretion forms around the compact object. The 
accretion disk prevents the outflowing material during the active phase and outflow with two-component velocity 
regime forms. 

2. The systems with parameters close to Z And have similar accretion rate of their compact object. The accretion 
rate determines the regime of hydrogen burning. The view that hydrogen burns in a steady state in the classical 
symbiotic stars is commonly accepted. When the accretion rate goes above the upper limit of the steady burning 
range the white dwarf expands which is observed as optical outburst with typical duration of one year. The 
first outburst opening the active phase can be followed by repeated outbursts. During the repeated outbursts an 
extended envelope surrounding the accretion disk exists in the system and it can be responsible for the collimation 
of the stellar wind. 



3. OBSERVATIONAL EVIDENCES IN SUPPORT OF THE MODEL 



3.1. The first outburst 



According to our model a geometrically thin accretion disk exists in the system in the beginning of the first outburst 
and the inner region of this disk can be optically thick. 

The presence of a moderate-velocity wind is indicated by the P Cyg profiles of the triplet He I A 4471 and He I 
A4713 lines during the light maximum. The two-component structure of their absorption showed mass outflow from 
the compact component with a mean velocity of 60 km s _1 (Fig. 3) [8]. 
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FIGURE 3. The profiles of the He I triplet lines in quiescence (January 1999; dot-dashed line) and near the maximum of the first 
outburst (November 2000; solid line). (From spectroscopy presented in Tomov et al. [8].) 

During the time of the increased light the Hy line had two emission components — a central narrow one and a broad, 
low-intensity component [8]. In November 2000 an absorption feature was probably present to the short-wavelength 
side of the narrow component, forming a dip in its emission profile (Fig. 4 right panel). The position of the dip 
corresponded to the red absorption component of the P Cyg lines of He I (Fig. 5). The low velocity of the P Cyg 
absorption and the dip in the H r emission profile can be explained with supposition for decrease of the velocity of the 
wind of the compact component after its collision with the disk (Fig. 1(b)). 
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FIGURE 4. The profile of the Hy line in quiescence of Z And, on January 7, 1999 (left) and at various stages of the 2000 
outburst (right). The dot-dashed line in the right panel corresponds to the spectrum on November 17, 2000, and the solid line — to 
the spectrum on December 5, 2000. The Gaussian fit of the broad component is marked with a thick line. The level of the local 
continuum is marked with a dashed line. (From spectroscopy presented in Tomov et al. [8].) 
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FIGURE 5. The profiles of the H y , He I A4471, and He I A4713 lines obtained on November 17, 2000. The positions of the two 
components of the He I absorption are marked with vertical dotted lines. (From spectroscopy presented in Tomov et al. [8].) 



The blue wing of the broad Hy component was not seen in November 2000 (Fig. 4, right panel). Assuming that the 
broad component (like the broad component of the He II A4686 line) is emitted by an optically thin high-velocity wind 
from the compact object, this can be explained with absorption in the stellar wind which is responsible for the P Cyg 
profiles. 



In the quiescent state an additional emission component was present on the short-wavelenghts side of the H r line 
forming a shoulder in its profile. Thus the full width at the zero intensity (FWZI) of the line reached to about 500 
km s _1 (Fig. 4, left panel). The presence of such a component can be explained with emission by an accretion disc. 

— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 




-1000 -500 500 1000 

Heliocentric RV (km/s) 



FIGURE 6. The wings of the He II 14686 line during the 2000 outburst (November 17). The broad component was visible in this 
period. (From spectroscopy presented in Tomov et al. [8].) 

The He II 14686 line of Z And was two-component one consisting of a central narrow component and a broad 
component with a low intensity [8]. The broad component was very weak in quiescence and could not be measured 
accurately. It became stronger during the first outburst and towards the light maximum its FWZI exceeded 1000 km s~' 
(Fig. 6). We believe that during this period the broad component of the He II 14686 and H r lines emitted in a region 
of stellar wind with a velocity of 500 kms -1 traveling at higher stellar latitudes. 

3.2. The fourth outburst 

The fourth outburst was the strongest one. The spectral data acquired during this outburst contain indications of all 
kinds of stellar wind, observed during the first one, and of bipolar collimated outflow in addition. We suppose that a 
disc-like envelope surrounding the accretion disc and collimating the outflowing material existed in the system at that 
time and the mass-loss rate was high enough to give rise to satellite line components. 
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FIGURE 7. The broad component of the H r (left panel) and Hell 14686 (right panel) lines during the outbursts in 2000 
(December 5; the dot-dashed line) and 2006 (October 4; the solid line). The level of the local continuum is marked with a dashed 
line. 

The analysis of our observations of October 2006 with high resolution (0.2 A per pixel, as for the first and second 
outbursts) in the H r and He II 14686 region of the spectrum of Z And reveals the presence of emission up to 500 
kms -1 from the centers of these lines (Fig. 7). Skopal et al. [9] obtained the effective photospheric radius of the 



outbursting component of (12 ± 4)Rq at the time of the light maximum for a distance of 1 .5 kpc. The minimal radius 
for the distance of 1.12 kpc used by us [8, 10] is 6Rq. If the broad emission originated in a Keplerian accretion disk 
at a distance equal to the radius of the outbursting component, its velocity should be not greater than 140 kms -1 , in 
disagreement with the observations. Like in the case of the first outburst it is natural to suppose that the observed broad 
emission is formed in the high-velocity wind at higher stellar latitudes (Fig. 2). A dip is observed in the blue part of 
the broad H r emission at a velocity of about 100 km s _1 . 
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FIGURE 8. Evolution of the H a line during the 2006 outburst. (From spectroscopy presented in Tomov et al. [12].) 

In our spectra [12], the H a line possessed a strong central component, a broad component (wings) and additional 
absorption and emission features to either side of the central component (Fig. 8). The central component was single 
peaked having a shoulder on its short-wavelength side, which was not seen only in the spectra taken in September. 
A small peak shortward from the central component was seen in the spectrum obtained on August 8, 2006, and the 
dip before the peak had a velocity about 100 kms -1 , in agreement with the results of the analysis of the Hy and He I 
profiles. Satellite emission components on either side of the central peak appeared after mid-July (Fig. 8). They are 
clear signs of the presence of a bipolar collimated wind. According to our model the emitting region is near the cone 
surfaces (see Fig. 2), where the density of the emerging collimated outflow is the highest. 

The first spectra during this outburst taken in July 2006 (Figs. 8 and 9) show a sharp absorption line at a velocity 
of 1400 kms -1 on the short-wavelength side of the central component and only a weak irregularly shaped emission 
component corresponding to a velocity of about 1500 km s _1 on its long-wavelength side. The blue-shifted absorption 
probably indicates the flow which is projected onto the stellar disk. The evolution of the spectrum in Figs. 8 and 9 
shows that the blue-shifted absorption component disappears and an emission appears. Two emission components on 
either side of the central peak were formed in July. The disappearance of the blue-shifted absorption component and 
the appearance of emission are probably due to a decrease of the mass-loss rate and/or an increase of the number of 
emitting atoms in the region of the wind that does not project onto the stellar disk. 
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FIGURE 9. The variation of the H a line: transition from absorption feature to emission one. The solid line corresponds to the 
spectrum taken on July 9, 2006; the dot-dashed line to July 14, 2006; the dotted line to August 8, 2006; and the dashed line to 
September 7, 2006. The level of the local continuum is marked with a dashed line. (From spectroscopy presented in Tomov et al. 
[12].) 



During all the outbursts considered by us the H a line possessed broad wings extending to velocities of at least 
±2000 km s -1 from its center [12, 8, 10]. Skopal et al. [7] concluded that the H a wings during the outburst state were 
formed in a high-velocity stellar wind from the compact component. However, as it was noted in Tomov et al. [10] 
it is not ruled out that the high-velocity broadening of the H a line was due to radiative damping and the stellar wind 
contributed to the emission at lower velocities. 



4. MASS-LOSS RATE 

The mass-loss rate of the compact object during the existence of the collimated wind was calculated from the energy 
flux of the broad Hy component, observed on October 3 1 , 2006 which was analyzed by fitting with a Gaussian function 
(Fig. 10). This broad component was not symmetric since its blue wing was absorbed by the wind outflow responsible 
for the P Cyg absorption components like during the first outburst. Having in mind this absorption we consider the 
wind velocity obtained by us as one lower limit. The equivalent width of the broad component was obtained with an 
error of 10 per cent depending mainly on the level of the local continuum. The real error however, is greater since the 
broad emission is not clearly visible because of blending with the central narrow component and its profile is actually 
not Gaussian. The energy flux of the broad component was obtained using its equivalent width and the continuum flux 
at the position of the line Hy. This continuum flux was calculated with linear extrapolation of the B and V fluxes from 
the photometric data of [20] taken very close to the time of our observations. The uncertainty of the continuum flux is 
not more than 10 per cent. The line flux was corrected for an interstellar extinction of E(B — V) = 0.30 according to 
the approach used in Cardelli et al. [21]. 

The mass-loss rate was calculated supposing that the medium of the wind is optically thin and that the outflow 
has a spherical symmetry and a constant velocity. The particle density in the wind is expressed via the continuity 
equation. In our calculations, we adopted a value of the electron temperature in the wind of 30000 K like during the 
first outburst [8]. We used a parameter fi — 1.4 [22] determining the mean molecular weight fim^ in the wind and a 
helium abundance of 0.1 [23]. We adopted a distance to the system d — 1.12 kpc [16, 5] as it was used in the paper of 
Tomov et al. [8] to compare the results more easily. It is supposed that the line is emitted by a spherical layer and the 
radii of integration must be estimated. We assumed optically thin medium and the inner radius in this case is thought 
to be the photospheric radius. The photospheric radius was estimated from the bolometric luminosity and the effective 
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FIGURE 10. The broad component of the H r line on October 31, 2006. The level of the local continuum is marked with a dashed 
line. 



temperature of the outbursting compact object at the time of our observation. We used a bolometric luminosity of 
10 4 Lq [9] and Zanstra temperature of 90000 K [13] and obtained a photospheric radius of 0.31(t//1.12kpc)/?Q. The 
outer radius of integration is equal to infinity. We used a recombination coefficients for case B [24] corresponding to 
temperature of 30 000 K and the density at the level of the photosphere. Having a line flux of 7.83 x 10~ 12 ergcm~ 2 s _1 
and a velocity of the wind of 370 kms -1 we obtained a mass-loss rate of 1.33 x 10~ 7 (<i/1.12 kpc) 3 / 2 MQ yr _1 . 

The total mass of the two parts (jets) of the collimated wind where the emitting and the non emitting masses are 
included, is equal to the mass, ejected by the system during the time of existence of the jets (July-December 2006). 
The ejected mass is m = MAt , where M is the mass-loss rate of the outbursting compact object and At — the time of 
existence of the jets. Then the total mass of the two jets is obtained to be about 0.7 x 10~ 7 (<f/1.12 kpc) 3 / 2 MQ. 



5. CONCLUSIONS 

The comparison of the results of our observations with the results of the gas-dynamical modeling provides possibility 
to propose a model to explain the behavior of the Z And system during its whole 2000-2010 active phase as well as 
the differences in its observed properties during different outbursts. The basic points of this model are the following 
ones: 

• An accretion disk exists in the system. 

• During the outburst the gas leaving the surface of the compact component at high velocity collides with the 
accretion disk. As a result its velocity decreases near the orbital plane, but remains unchanged at higher latitudes. 
Hence, two wind components with different velocities are observed. 

• During the outburst the wind from the compact component partially disrupts the disk. Some portion of its mass 
does not leave the potential well, after the wind ceases it begins to accrete again. Because of the initial amount of 
angular momentum a disk-like envelope extending to larger distance from the orbital plane forms. 

• The existence of a centrifugal barrier creates hollow cones around the axis of rotation which lead to appearance 
of collimated outflow during the recurrent outbursts. This outflow can be observed as high-velocity components 
on either side of the main peak of the line. 

• Collimated outflow will appear when the density of the disk-like envelope is high enough to provide collimation 
and the mass-loss rate of the compact component is also high. High-velocity line components thus can be observed 
during outbursts when the loss of mass is at high rate, preceded by a similarly strong outburst. 

The model suggested provides possibility to explain all the spectral features of Z And observed during active stage. 
The general character of this model permits us to suppose that similar scenario is possible for other classical symbiotic 
stars too. 
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